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Abstract 



We study whether spin-independent scattering of weakly-interacting massive par- 
ticles (WIMPs) with nuclei can account for the annual modulation signal reported by 
DAMA. We consider both elastic and inelastic scattering processes. We find that there 
is a region of WIMP parameter space which can simultaneously accommodate DAMA 
and the null results of CDMS, CRESST, and XENON. This region corresponds to 
an ordinary, elastically-scattering WIMP with a standard Maxwell-Boltzmann distri- 
bution, a mass 3 GeV < uidm i$ 8 GeV, and a spin-independent cross section with 
nucleons 3 x 10~ 41 cm 2 < er^ < 5 x 10~ 39 cm 2 . This new region of parameter space 
depends crucially on the recently discovered effect of channeling on the energy thresh- 
old for WIMP detection in the DAMA experiment; without the inclusion of this effect, 
the DAMA allowed region is essentially closed by null experiments. Such low-mass 
WIMPs arise in many theories of Beyond the Standard Model physics, from minimal 
extensions of the MSSM to solutions of the baryon-dark matter coincidence problem. 
We find that inelastic scattering channels do not open up a significant parameter region 
consistent with all experimental results. Future experiments with low energy thresholds 
for detecting nuclear recoils, such as CDMSII-Si and those utilizing ultra-low energy 
germanium detectors, will be able to probe the DAMA region of parameter space. 



1 Introduction 



Recently, the DAMA collaboration has provided further evidence for the observation of an 
annual modulation in the rates of nuclear recoil in their experiment [1]. Such a signal arises 
naturally from postulating Weakly Interacting Massive Particles (WIMPs) in the galactic 
halo that scatter from target nuclei in detectors. The annual modulation of the interaction 
rate comes from the variation in the relative velocity of the earth with respect to the galactic 
dark matter halo as the earth orbits the sun. This changes the flux of dark matter particles 
and the size of their interaction cross-sections, with expected extrema occurring at June 2 and 
December 2. The DAMA experiment observes a maximum at low nuclear recoil energies on 
May 24, plus or minus 8 days, and they have accumulated enough data to put the significance 
of the observed modulation at approximately 8a. Both the phase and amplitude of the signal 
are highly suggestive of WIMP interactions. The collaboration has not been able to identify 
other systematic effects capable of producing this signal, and have claimed that the annual 
modulation is a discovery of dark matter. This claim has been controversial, partly because 
a number of other experiments seem to be in direct contradiction. In particular, the original 
DAMA allowed region with WIMP mass 30 GeV < m DM < 200 GeV and dark matter- 
nucleon interaction cross-section o v ~ 10~ 41 — 10~ 42 cm 2 had been quite conclusively ruled 
out by the CDMS [2-4] and XENON [5] experiments for the case of an elastically scattering 
WIMP. 

Methods of reconciling the DAMA signal with the results of other experiments have 
been proposed in the past. Inelastic scattering processes, Xi-^ — * X2N, where Xi is the 
dark matter particle and \2 is another new state with mass splitting 5 between the states, 
have been proposed in the context of supersymmetric models [6-8]. Inelastic scattering 
of MeV dark matter particles to lighter states was investigated as a possible solution in 
Ref . [9] . Mirror states from a hidden-sector copy of the Standard Model have been proposed 
as a candidate consistent with all experimental constraints [10, 11], as have various models 
with heavy composite states [12]. A model-independent study of spin-independent elastic 
scattering noted that scattering from the sodium component of the Nal DAMA scintillators 
allowed a small window of dark matter masses in the 5 — 9 GeV to be consistent with 
current experimental constraints [13, 14]. This study noted that since sodium nuclei are 
lighter than germanium nuclei, the threshold for scattering off sodium could be lower than 
that for germanium for light dark matter states. It was also shown previously that spin- 
dependent scattering may open up additional parameter space consistent with DAMA and 
other experiments [15]. 

This mo del- independent study of elastic scattering in Ref. [14] is no longer applicable, 
as several new results have appeared recently. New experimental constraints and improved 
understanding of scattering processes in the DAMA apparatus have drastically altered both 
the excluded parameter space and the physics underlying the DAMA modulation signal. 
We perform a model-independent study of both elastic and inelastic scattering mechanisms 
accounting for all recent experimental measurements. We find that completely ordinary, 
spin-independent elastically scattering WIMPs with masses in the range 3 — 8 GeV and 
scattering cross sections in the range 3 x 10 -41 cm 2 to 5 x 10~ 39 cm 2 are consistent with all 
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experimental constraints. No additional dark matter stream is needed; a simple Maxwell- 
Boltzmann distribution allows this parameter space. Inelastic scattering no longer opens up 
a significant region of additional allowed parameter space. We summarize here the important 
features and conclusions of our analysis. 

• We include the effect of channeling in the Nal crystal scintillators of DAMA, an effect 
recently noted in Ref. [16] and studied by the DAMA collaboration [17]. Channeling 
occurs in crystalline detectors where the only signal measured is the light output, and 
when recoiling nuclei interact only electromagnetically with the detector material be- 
cause of either their direction of motion or incident energy. The effect of channeling 
is to remove the quenching factor usually required to convert between nuclear recoil 
energy and electron-equivalent energy, and in the context of dark matter searches it 
effectively lowers the energy threshold for detection of nuclear recoils of DAMA below 
that of CDMS and XENON. This effect is crucial in reconciling elastically scattered 
WIMPs consistent with all experimental constraints. In particular, the lower thresh- 
old of DAMA means that it can detect lighter dark matter particles than the higher 
threshold experiments like CDMS and XENON. This effect opens the region of light 
WIMP parameter space for DAMA. We note that the presence of channeling in the 
energy regime studied by DAMA has not been conclusively established, although it 
has been observed in Nal crystals at higher energies [18]. 

• We include constraints from CDMS-SUF, CDMS-II, CRESST-I scattering from sap- 
phire targets [20] , and XENON. Inclusion of experimental results from multiple target 
nuclei is necessary to correctly elucidate the allowed parameter region. In the elastic 
scattering case we also include the recent results from the CoGeNT collaboration [19]. 

• We study spin-independent elastic scattering, and inelastic scattering with either posi- 
tive or negative mass splitting 6 between the incident and scattered dark matter parti- 
cle. Inelastic scattering of either sign opens up only a very small window of parameter 
space; roughly, inelastic scattering to heavy states is ruled out by XENON and the ger- 
manium data from CDMS-II, while scattering to lighter states is ruled out by CRESST 
and the silicon data from CDMS. The preferred parameter space is for light mass dark 
matter with elastic scattering. Future results from ultra-low noise germanium detec- 
tors [21] 0, and the lower threshold silicon data from CDMS will be vital in exploring 
this region. The low threshold germanium experiment TEXONO [22] may also be able 
to probe this region, though its sensitivity has recently been called into question [23] 
(inclusion of current TEXONO constraints does not change our results). 

There are many possible models which could give rise to such a comparatively light 
WIMP. In extensions of the Minimal Supersymmetric Model, for example, GeV mass WIMPs 
with the right relic abundance arise [24,25]. It has been shown that hidden sectors in the 
context of supersymmetric models give rise naturally to WIMPs with GeV or even lighter 

1 We thank J. Collar for correspondence regarding the ability of these experiments to probe this region. 
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masses, as observed in Refs. [26,27]. It was shown explicitly that the models of this sort 
can account for the DAMA signal [28]. Supersymmetric models with non-unified gaugino 
masses at the grand unified scale give rise to light neutralino dark matter candidates [29]; 
the importance of the channeling effect for these models was noted in Ref. [30], as was 
the effect of having the light sodium component of the DAMA target. Lastly, solutions 
to the baryon-dark matter asymmetry problem also predict a WIMP with a mass in the 
range, m^ m ~ f2d m /f2&m p ~ 5m p [31-34]. It is clear that the light dark matter paradigm 
suggested by the direct detection experiments raises numerous theoretical questions and 
has phenomenological impact on a broad array of experiments. We leave the potential 
implications of light WIMPs to future work, and focus here on clarifying the experimental 
situation. 

The outline of the paper is as follows. In Section [2] we review the formalism of direct 
detection of dark matter. We review the characteristics of the relevant experiments, discuss 
the physics and implications of the channeling effect in DAMA, and discuss our analysis 
method. We then apply these techniques to derive the allowed parameter space consistent 
with all experimental measurements for both elastically and inelastically scattered WIMPs 
in Section [31 Finally, we conclude and discuss future directions. 



2 Rates for direct detection of dark matter 
2.1 Direct detection formalism 

We first review the relevant features of direct detection of dark matter by WIMP-nucleus 
scattering. The differential rate per unit detector mass in nuclear recoil energy is given by 

dR Ar Pdm f .;; ,., , ., da 



ir jVy— — — — / d*vvf(v,v e )^-. (2.1) 

Here, Nt is the number of target nuclei per unit mass, ttidm is the dark matter particle mass, 
and pdm = -3 GeV/cm 3 is the local dark matter halo density, v is the dark matter velocity 
in the frame of the Earth, v e is the velocity of the Earth with resepect to the galactic halo, 
and f{v,v e ) is the distribution function of dark matter particle velocities. We take f{v,v e ) 
to be a standard Maxwell-Boltzmann distribution: 

MO = ^4^ e ~ (tmre)> ° 2 - ( 2 - 2 ) 



[rev, 



2)3/2 



The Earth's speed relative to the galactic halo is v e = v Q + t> or f>cos7 cos[u(t — t )] with 
v q = v o + 12km/s, v or b = 30km/s, cos7 = 0.51, t Q = June 2nd, and uj = 27r/year. We set 
the most probable dark matter speed in the galactic frame to vq = 220km/s for most of our 
analysis, but study the effect of allowing it to vary within its 90% C.L. range 170km/s < 
vo < 270km/s [35]. The upper limit of the velocity integration of Eq. ( 12. ip should be taken 
as the galactic escape velocity, 490km/s < v esc < 730km/s at 90% C.L. [35]. For simplicity 
we take the upper limit to infinity and a postieri restrict the minimum velocity to lie above 
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v esc . We set v esc = 730km/s for our study, but again study the effect of varying it within 
the 90% CL. range. Neither uncertainty has a significant effect on our conclusions. The 
minimum dark matter velocity v m i n depends on the recoil energy E R and the details of the 
direct detection process, and will be discussed later. 

We assume a spin-independent cross section between dark matter particles and nuclei; a 
standard calculation [36] leads to 

da _ m N a p [f p Z + f n (A - Z)] 2 p2 ^ 



dE R 2v 2 /4 ft 

fi n is the reduced mass of the dark matter particle and nucleon (proton or neutron), a p 
is the scattering cross section of the dark matter particle with nucleons, and f UtP are the 
coupling strengths of the dark matter particle to neutrons and protons respectively. f n<p 
are calculated from a coherent sum over the couplings to the quark model constituents of 
the nucleon. Z and A are the proton and atomic numbers of the nucleus, while F(E R ) is a 
nuclear form factor which corrects for the above formula being strictly correct only as the 
momentum transfer q 2 = 27jinEr — ► 0. We use a standard Helm form factor [37]. 

Inserting the dark matter-nucleus scattering cross section of Eq. (12.31) into the differential 
rate in Eq. (12.11) and performing the integration, we obtain the differential rate 

dR N T m N p DM a p [f p Z + f n (A - Z) 2 } 2 f erf (x min + rj)- erf (x min - rj) 



dE R Av m DM /4 fn 

x9(v e sc-v miri (E R )) . (2.4) 

We have introduced the parameters x min = v min /vo and 77 = v e /v . The 6 function accounts 
for our treatment of the maximum allowed dark matter velocity, as discussed above. We will 
use the formula in Eq. (12.41) to interpret the results of various direct detection searches for 
dark matter particles. We follow the standard convention and set f p = f n ; results for other 
exchange mechanisms, such as Z-boson exchange, can be derived by scaling the cross section 
a p appropriately. The total number of nuclear recoil events in a recoil energy range between 
El and E 2 is 

N = Y r Er dR^£i dE (2 5) 

^W £ i dE r N T Mi 1 ; 

where the sum is over each nuclear species i in the detector. Si is the effective exposure 
of species i expressed in kg-days. We can parametrize Si = Mrfe, where t is the time of 
exposure, Mj is the target mass of species i, and e is a detection efficiency. 



2.2 Experimental constraints 

We discuss the direct detection signal obtained by DAMA and the constraints imposed by 
other null experiments. The DAMA collaboration has recently reported an 8.2a signal for 
an annual modulation signature of dark matter scattering [1] . This signal has been reported 
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for several years, and the significance of the effect has been constantly increasing. It is 
convenient to parametrize the DAMA detection rate as 

Ri = 1$ + S}cos[u(t -t )], (2.6) 

where to = June 2nd = 152.5 days is the time when the Earth is moving with its maximum 
speed with respect to the galactic halo and T = 2ti/u = lyr. The subscript i denotes 
different energy bins. The values measured by the collaboration in the 2 — 6 keVee bin are 
to = 144 ± 8 days and T = 0.998 ± 0.003 yr. The keVee unit is keV-electron-equivalent, 
which accounts for quenching of unchanneled events; we describe this effect in detail in the 
next subsection. For the moment, we only comment that keVee = keV for channeled events, 
and keVee = qkeV, where q < 1, for unchanneled events. The constant term R® is composed 
of both a signal piece coming from dark matter initiated processes, and a background piece 
arising from other sources of nuclear recoil: = + Sf. The expressions for S® and Si 
are obtained by integrating Eq. (12.41) over a given range of Er. More precisely, expanding 
the interaction rate in a Taylor series around v e = v & for an energy interval between Ei and 
Ei + AEi, we have 

and 

_ Av e f^' a (dR\ 

where the differential rates are to be evaluated at v e = v Q and Av e = v or bCos , y. 

The DAMA collaboration has not been able to find another effect besides dark matter 
scattering that could contribute to S}. The modulation amplitudes reported by DAMA in 
each energy bin are given in Table [TJ Non-zero modulation is clearly observed in the low 
energy bins, while it is absent in the 6 — 14 keVee bin. 



Energy 


Sj (cpd/kg/keVee) 


2-4 keVee 


0.0223 ±0.0027 


2-5 keVee 


0.0178 ±0.0020 


2-6 keVee 


0.0131 ±0.0016 


6-14 keVee 


0.0009 ±0.0011 



Table 1: Modulation amplitudes in units of counts-per-day/kilograms/keVee as reported by 
DAMA [1] 

Any dark matter model allowed by the DAMA results must be consistent with other 
null experiments. We consider constraints imposed by the CDMS, CRESST, and XENON 
collaborations. We summarize the salient features of each experiment, including the target 
nucleus, energy threshold, and exposure, in Table [2j Considering the constraints from a 
variety of experiments is vital, as each is relevant for different regions of dark matter pa- 
rameter space. We will see this explicitly in our numerical results of Section [3j The most 
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important qualitative observation to make about the experiments is the difference in the 
energy thresholds. The recoil energy is related to the mass of the incoming dark matter 
particle through the relation 

E R = ^Jfaac, (2.9) 
m N 

where fi nuc the reduced mass of the WIMP- nucleus system and v their relative velocity. It 
follows that for low mass WIMPs, the recoil energy of the nucleus is not high enough to 
be above threshold for XENON and CDMS, while it is for DAMA. The experiments which 
do have thresholds competitive with DAMA, in particular the CRESST-I experiment, have 
much lower exposures. While they constrain the low dark matter mass region, their limits 
on the cross-section a p are not strong enough to constrain much of the DAMA parameter 
space. We show in Section [3] that the CDMS-II Si run is also sensitive to lower dark mater 
masses because of the lower target mass of silicon and the lower threshold compared to the 
germanium data. 



Experiment 


Target 


Exposure (kg-d) 


Threshold 


Ref 


CDMS-SUF 


Ge 


65.8 


5 keV 


[2] 




Si 


6.58 


5 keV 




CDMS-II 


Ge 


121.3 


10 keV 


[3] 




Si 


12.1 


7 keV 


[4] 


XENON10 


Xe 


131 


4.5 keV 


[5] 


CRESST-I 


A1 2 3 


1.51 


0.6 keV 


[20] 



Table 2: Relevant features of the null experiments used in our analysis. 



2.3 The physics of quenching and the DAMA signal 

Since it has important implications for dark matter searches, we discuss in more detail the 
physics of the channeling and quenching effects which allow the compatibility of the DAMA 
signal with other null experiments. The relevance of this effect for the current analysis is 
that DAMA is sensitive to nuclear recoil energies down to 2 keV, below the threshold of both 
XENON and CDMS. 

The light yield of scintillators is different depending on whether the incident particles 
interact electromagnetically or via the nuclear force. The bookkeeping for this is performed 
by introducing an electron equivalent recoil energy for a given nuclear recoil energy; for 
a nuclear recoil of 1 keV, the equivalent electron energy (in keVee) is q x x 1 keV. q x is 
the quenching factor for the nuclear material composing the scintillator; for the materials 
composing the DAMA detectors, q^a ~ 0.3 and qj ~ 0.09. It is simple to understand 
why q < 1 and the electron equivalent energy is less than that of the nuclear recoil energy. 
A nucleus hitting scintillator material will lose energy both by electromagnetic and nuclear 
interactions, while an incoming electron will lose energy only via electromagnetic interactions. 
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It is primarily the production of radiation in electromagnetic interactions that produces the 
light yield in the scintillating material, yielding q < 1. 

For crystal scintillators such as those used by DAMA, a portion of the events will be 
"channeled," effectively changing the quenching factor to q ~ 1. This occurs when incident 
particles that do possess nuclear interactions interact only electromagnetically with the scin- 
tillator material. This can occur for certain energies and incidence angles of the incoming 
particle. The importance of this for the DAMA experiment was first discussed in Ref. [16], 
and a detailed analysis of its effect was performed by the DAMA collaboration in Ref. [17]; 
we refer the reader to these references for a more detailed discussion. An estimate of the 
fraction of channeled events based on simulation is given in Fig. 4 of Ref. [17]. We use the 
following simple parametrization in our analysis [11]: 

fNa ~ l + l.UE R (keVY fl ~ 1 + 0.75E R (keV) ' (2 ' 10) 

The other experiments considered in our study do not report this distinction between 
quenched and channeled events. CDMS collects all recoil energy from ionization and phonons 
(heat). Channeling does not occur in liquid noble gases [17], such as used in XENON. 
CRESST sets q — 1 in their analysis. The XENON collaboration has incorporated the 
appropriate quenching factor into their published results. 



2.4 Analysis procedure 

We now discuss our analysis procedure for determining whether the DAMA results are con- 
sistent with the constraints of other experiments for both elastic and inelastic dark matter 
scattering. For experiments that report no events above background, we demand that dark 
matter initiated scattering produce less than 2.3 events throughout the entire range analyzed, 
as appropriate for establishing a 90% C.L. upper limit. For experiments that report events 
but ascribe them to background processes, we adopt a conservative approach and include 
the events when setting limits. We determine the 90% C.L. limit on the dark matter cross 
section in the presence of this background using a simplified version of the optimum interval 
method [38]. We divide the data into energy bins using the published results. We study all 
contiguous combinations of bins, and demand that the dark matter candidate produce no 
more than the 90% C.L. upper limit allowed events as dictated by Poisson statistics in each 
such interval. The most stringent limit in a p for a given uidm among the studied intervals is 
then used. This procedure allows effective use of the different kinematics of each dark mat- 
ter particle when setting constraints, as the appropriate energy bins which constrain a given 
dark matter candidate depend strongly on the mass and whether elastic or inelastic scatter- 
ing is considered. Our technique reproduces reasonably well the published constraints in the 
elastic scattering limit, as obtained from either the experimental papers or DM Tools [39]. 
For completeness, we now describe in detail our treatment of each experimental data set. 

• DAMA : For the DAMA signal, we perform a \ 2 fit to the 2 — 6 and 6 — 14 keVee bins 
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given in Tabled] for a given ttidm and cr p : 




2 



(2.11) 



where R\ ,exp is the value of R\ measured by DAMA in the ith energy bin, R\ is the 
rate computed at a given (m DM , a p ) point, and <jj is the experimental error on R\ ,exp . 
The minimum x 2 is found at a given tudm by scanning over a p . Once the minimum 
X 2 is found, if Xmin < 2, a 90% C.L. region is determined by accepting values of o v 
with x 2 — Xmin < 2.71. We also determine a 3a allowed range for a p by accepting cross 
sections leading to x 2 ~ Xmin < 9- 

• CDMS-II : The latest results from the five-tower germanium run of CDMS indicate no 
events in a nuclear recoil energy range of 10 — 100 keV [3]. The efficiency is fairly flat, 
and it has already been included in the effective exposure given in Tabled We demand 
that a dark matter candidate produce < 2.3 events to establish a 90% C.L. upper limit 
on cr p for a given tudm- We also include the constraint from the silicon run of Ref . [4] . 
The energy window for this result is 7 — 100 keV; the lighter mass of silicon and the 
lower energy threshold leads to important constraints for light dark matter. We again 
demand that dark matter interactions produce < 2.3 events. 

• CDMS-SUF : We include constraints from the CDMS run reported in Ref. [2]. This 
early run has implications for light dark matter because the analysis threshold was set 
at 5 keV. The efficiency is not flat as a function of energy. We use the parametrized 
form given in Ref. [14]: 



This function must be included in the integration of Eq. (I2.5p . The experiment does 
observe events in the expected dark matter signal range, although these are consistent 
with background expectations. We determine the 90% C.L. limit on the dark matter 
scattering rate in the presence of this background using a simplified version of the 
optimum interval method discussed above. We divide the data into the four energy 
bins E < lOkeV, lOkeV < E < 20keV, 20keV < E < 55keV, and E > 55keV. We 
obtain the background events in each bin from Ref. [2]. 

• XENON : We include the first results from the XENON dark matter search reported 
in Ref. [5]. We again apply the optimum interval method using the efficiencies and 
energy bins given there, and the reported background events. All observed events are 
attributed to background processes by the collaboration, but we follow their analysis 
in including the observed events when computing a limit. 




E < 10 keV 
10 keV < E < 20 keV 
E > 20 keV. 



(2.12) 
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• CRESST : The initial run of the CRESST experiment utilized a sapphire target with 
a low threshold of 0.6 keV [20]. We use the observed spectrum without background 
subtraction given in that reference. The collaboration attributes all observed events 
to background initiated processes, due to a study of coincident counts in multiple 
detectors. We apply the optimum interval method with a restriction that the interval 
size be greater than 1.2 keV, which is twice the experimental energy resolution. We 
use a 100% efficiency for detection of events, as reported by the collaboration. 

• TEXONO : The TEXONO collaboration has recently reported results using a germa- 
nium target with a 200 eVee nuclear recoil energy threshold [22] . The limits claimed in 
this paper (but questioned in [23]) do not yield significant improvement over CRESST- 
I. We have checked that the 3 — 8 GeV WIMP mass window is still open with or without 
the inclusion of TEXONO constraints in the current analysis. Their projected sensi- 
tivity, however, cuts deeply into the DAMA parameter space. We include the future 
projections in our analysis. 

• CoGeNT : In the elastic scattering case we include the recent result of the CoGeNT 
collaboration [19]. The current exclusion curve reported by the collaboration does not 
include the possible effect of channeling. It is currently unclear by how much Ge- 
based ionization detectors such as used in CoGeNT are affected by channeling. We 
caution the reader that the CoGeNT exclusion curves could drift on our plots when 
channeling in this type of detector is simulated. Without channeling factored in, the 
measurement does not close the light WIMP window opened up by the channeling 
effect in the DAMA experiment, although future results from the collaboration should 
either confirm or severely constrain the DAMA parameter region. 

3 Implications for models of WIMP dark matter 

Following the analysis procedure described in the previous section, we study whether pa- 
rameter space exists in which spin-independent dark matter scattering off nuclei can si- 
multaneously accommodate the DAMA modulation signal and satisfy the constraints from 
other experiments. We begin by considering the case of elastic scattering. We then examine 
whether inelastic scattering processes of the form \iN — > X2N, where \2 can be either lighter 
or heavier than the dark matter particle xi> can account for all experimental observations. 

3.1 Elastic scattering 

We present in Fig. ([1]) the results of our parameter space scan for spin-independent dark 
matter scattering. The inner and outer hatched regions respectively denote the 90% C.L. 
and 3a regions consistent with the DAMA modulation signal. Regions above each colored 
line are excluded by the indicated experiments. Although a portion of the region consistent 
with DAMA is excluded, a region with w,dm ~ 3 — 8 GeV is not excluded. Two features 
of current experiments make this possible. Substituting the 4.5 keV threshold for XENON 
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and the 10 keV threshold for CDMS-II into Eq. (12. 9p and setting v = v esc = 730km/s, we 
see that XENON can only probe tuom > 7GeV, while CDMS-II germanium is sensitive to 
TfiuM > 8 GeV. Although these experiments have large exposures, the light dark matter 
region is inaccessible to them. Further running of these experiments will not completely 
remove this window. The experiments with low enough thresholds to be sensitive to the 
DAMA region, notably CRESST, the 7 keV CDMS-II silicon run, and CoGeNT, do not have 
a large enough exposure, as seen from Table [2j 

Several distinct structures are noticeable in the DAMA allowed region. The low mass 
tail extending to large cross sections arises from channeled sodium events. Both the low 
threshold and light nucleus cause these to appear first as uidm is increased. The allowed 
region above m^M ~ 50 GeV arises from quenched iodine events. The presence of channeling 
is crucial for reconciling the DAMA signal with other null experiments. This is illustrated in 
Fig. (J2]), where the effects of channeling have been removed. No region in parameter space 
is consistent with all experiments when channeling is not considered. 

The previous analysis of elastic scattering in Ref. [14] utilized the light mass of sodium 
to render DAMA consistent with other experiments. At the time, all events in DAMA were 
considered to be quenched, leading to an approximately 2 keV/qj ~ 7 keV energy threshold 
for scattering off sodium. Although this threshold is similar to the thresholds of other 
experiments, sodium is lighter than other nuclei such as germanium, so that the net effect 
is sensitivity to lighter dark matter particles. This led to the low-mass tail in Fig. (T5]). New 
limits from the 7 keV silicon run of CDMS-II have completely ruled out this region, although 
as shown in Ref. [14] it is possible that fitting to the 2 — 4 keV bin instead of the 2 — 6 keV 
bin might open up a small amount of additional parameter space. The channeling effect is 
required to open a significant low-mass elastic scattering window. 

To roughly study the effect of galactic uncertainties on these results, we also show in Fig. 
(12]) the effect of changing the most-probable speed v and the galactic escape velocity. In 
the left panel we change to vq = 170km/s, while in the right panel we set v esc = 610km/s. 
Although both changes have quantitative effects, the picture described above is unchanged. 
Previous analyses have shown that assuming a galactic velocity distribution beyond a simple 
Maxwell-Boltzmann form opens up a small region of additional allowed parameter space [14], 
and we expect the same would occur in our study, though we do not pursue it further here. 

Further running of these experiments will continue to probe the elastic scattering pa- 
rameter space. To estimate whether future experiments can rule out the allowed region, we 
show in Fig. (g]) the projected sensitivities of XENON, CDMSII-Si, and TEXONO together 
with the DAMA allowed region. These projections are obtained by simple scalings of the 
current exposure by a factor of ten. For CDMS-II, the expected signal has been increased by 
ten, while for XENON the expected signal and observed background have both been scaled 
upward. While the XENON energy threshold is currently too high to probe the DAMA 
allowed region, future CDMSII silicon data will test dark matter in the 5 — 8 GeV range. For 
TEXONO, we have assumed an exposure of 1 kg-yr, the future goal of the collaboration [22]. 
For CoGeNT, we have used the projected bound from Ref. [19] expected after an ongoing 
detector upgrade. Achievement of these goals will allow both CoGeNT and TEXONO to 
study the entire low-mass WIMP window. 



10 




12 5 10 20 50 100 



m DM (GeV) 

Figure 1: Allowed region in the mDMi&p plane consistent with the DAMA modulation 
signal at 90% CL. and 3a (inner and outer hatched regions, respectively). Also shown are 
the experimental constraints arising from other null experiments. The DAMA allowed region 
includes both channeled and quenched events. 




m DM (GeV) 

Figure 2: Similar to Fig. (JTJ), but if DAMA observed only quenched events. The presence of 
un-quenched (channeled) events is necessary to reconcile DAMA with null experiments. 
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Figure 3: Region allowed by the DAMA modulation signal together with other experimental 
constriants obtained by changing to vq = 270 k/s (left panel) and v esc = 610km/s (right 
panel). 



In summary, completely ordinary, elastically scattering WIMPs with masses in the 3 — 
8 GeV range can produce the DAMA modulation signature and remain consistent with other 
experimental constraints. Future data, particularly those from the silicon detectors of CDM- 
SII, from TEXONO and from other low-threshold Germanium experiments such as used by 
the CoGeNT collaboration, are critically important to probe this parameter region. 



3.2 Inelastic scattering 

Inelastic scattering has also been considered as a mechanism for explaining the DAMA signal. 
In models of this type, scattering processes xiN ~* X2N are considered. N denotes the target 
nucleus and %i is t ne dark matter candidate. The final-state particle \2 can have a different 
mass than xi- We denote the mass splitting by 5 = ni2 — mi. 5 is positive for heavier final 
states, such as in the supersymmetric model of Refs. [7,8]. Negative 5 has been studied as 
a possible phenomenological explanation of the DAMA signal in Ref . [9] . 

In the case of inelastic scattering, the maximum recoil energy can be shown by kinematics 
to be 

ET X = — l^iP 2 ~ 5 + Jii\P-2fr5p\ . (3.13) 
m N { v J 

Equivalently, for a given recoil energy, the dark matter must have a minimum velocity v m i n 
to scatter off a nucleus which is 

v mm = J^— (nt^ + s). (3.14) 

Here, fix = m\m^j{m\ + uin) is the reduced mass of the nucleus and Xi> an d (3 is the 
velocity of xi i n units of c. Inelastic dark matter with 5 > reconciled DAMA with 
the null experiments through two main effects [7,8]. First, for a sufficiently large mass 
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Figure 4: Projected sensitivities of XENON and CDMSII-Si obtained by scaling the current 
exposures upwards by a factor often. For XENON, both signal and observed background are 
scaled upwards, For TEXONO, the projected sensitivity has been read off from [22], while 
for CoGeNT the projected sensitivity is taken from [19]. 

splitting, only heavier targets can scatter inelastically for an experiment with a threshold 
E^. The restriction present in Eq. (1 3 . 1 3 j) that the dark matter velocity (3 2 > 25/ fix becomes 
increasingly stringent as fii is decreased. Since iodine is relatively heavy in comparison to 
germanium, this could enhance the signal in DAMA relative to the germanium of CDMS. 
This argument becomes much less effective with the presence of XENON constraints, since 
the mass of xenon is similar to iodine. Second, positive 5 changes the shape of the recoil 
energy spectrum, suppressing low energy recoils with respect to high energy ones. Since 
before the discovery of channeling it was thought that the thresholds of DAMA were higher 
than those of the null experiments, DAMA favored inelastic dark matter. 

However, we find that with the inclusion of channeled events and the new constraints 
from XENON and CDMS-II, inelastic dark matter becomes severely constrained. We show 
in Fig. <§5§ the DAMA allowed parameter space together with other experimental constraints 
for the representative values 5 = 25,50keV. The region of weak scale dark matter tudm ~ 
100 GeV and 100 keV mass splittings proposed in [7, 8] is closed. The new experimental 
constraints from XENON and the germanium run of CDMS-II are too severe. The small 
sliver of allowed parameter space present in the 5 = 25keV case arises from channeled 
scatterings from iodine. Although some regions for not-too-large 5 are still consistent with 
all measurements, inelastic scattering with 5 > does not open up significant parameter 
space beyond elastic scattering. 

There are two additional features to note in Fig. (J5J). First, relative to the 5 = case, 
positive 5 shifts the region allowed by DAMA to larger masses. This occurs because of the 
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higher incident energy required to push the scattering over threshold. Second, the constraints 
from experiments with light targets such as CRESST become much less significant. The 
constraint on the initial dark matter velocity f3 2 > 25 discussed before indicates that for 
light nuclei, only the high velocity tail of the Maxwell-Boltzmann distribution can lead to 
dark matter scattering. 




m DM (GeV) m DM (GeV) 

Figure 5: Region allowed by the DAMA modulation signal together with other experimental 
constriants in the m DM , a p plane for 5 = 25, 50keV. 

We now study the region 5 < 0. We display in Fig. (E]) the DAMA allowed pa- 
rameter space together with other experimental constraints for the representative values 
5 = — 25, — 50keV. The DAMA allowed region shifts to lower dark matter masses. The 
lighter-mass horseshoe-shaped region is associated with channeled events from sodium, while 
the region for slightly higher masses comes from both quenched sodium and channeled io- 
dine events. The lack of an allowed region at masses higher than roughly 10 GeV occurs 
because the phase of the DAMA signal cannot be correctly obtained; a minimum in the dark 
matter scattering rate is predicted at to — June 2nd rather than a maximum. Only a small 
sliver of the higher-mass region is consistent with all measurements, with the light-mass 
parameter space arising from sodium events ruled out by CRESST. The increasing severity 
of the these constraints as 5 is descreased can be understood by studying the regions of 
the Maxwell-Boltzmann distribution contributing to dark matter scattering as a function of 
5. Only a small portion of the Maxwell-Boltzmann distribution can initiate scattering for 
5 > 0. Lower velocities very quickly begin to lead to scattering as 5 is decreased, strength- 
ening the constraints for negative 5. The sharp turn-off of the CRESST constraints at low 
dark matter masses visible in Fig. (jH]) follows from the fact that Eq. (13.131) requires that 
m i <i m N E t ft Tesh J '(y lc 2m^E l ^ resh (3^ sc + for dark mater particles significantly lighter than 
the target nucleus. We note that the wiggles visible in the exclusion curves arise from our 
use of the binned optimum interval method; as the kinematics of the dark matter scattering 
changes, a different energy interval is selected to provide the most stringent constraint. 

Finally, we show in Fig. (171) the region in the mDM-,5 plane consistent with both the 
DAMA signal and all considered experimental constraints. Dark matter masses in the region 
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Figure 6: Region allowed by the DAMA modulation signal together with other experimental 
constriants in the tudm, Op plane for 5 = —25, — 50keV. 

3 — 13 GeV and mass splittings in the range — 15keV < 5 < 30keV are allowed. We note, 
however, that the precise width of these intervals is quite sensitive to small deviations in the 
constraint curves and in the signal. A combination of constraints from all null experiments is 
critical in obtaining this picture; XENON and CDMSII-Ge constraints are most important 
for positive 5, while CRESST and CDMSII-Si close off the negative 5 window. The largest 
range of allowed dark matter masses occurs for 5 ~ 0, indicating that inelastic scattering 
processes are not very helpful in reconciling the DAMA signal with other experimental 
constraints. 

4 Conclusions 

We have studied the consistency of the dark matter interpretation of the annual modulation 
signal observed by DAMA with the results of the null experiments CDMS, CRESST, and 
XENON. Recent work has shown the presence of a channeling effect in the crystal scintil- 
lators utilized by DAMA which drastically changes the interpretation of the experimental 
results. The presence of the channeling effect opens a window in dark matter parameter space 
between 3 and 8 GeV where the DAMA signal is consistent with all of the null experiments. 
This consistency requires no exotic dark matter physics-a vanilla, elastically scattering dark 
matter candidate interacting through spin-independent channels is sufficient to explain both 
the signal and the null results from the other experiments. We have also examined whether 
possible inelastic processes can accommodate all experimental results. Inelastic scattering of 
dark matter particles to heavier final states renders dark matter masses up to approximately 
13 GeV consistent with all measurements. However, the largest range of permissible dark 
matter masses occurs for elastic scattering candidates, indicating that inelastic processes 
do not open up significant regions of parameter space. Future measurements from ultra- 
low energy germanium detectors and silicon results from CDMS are needed to explore this 
light-mass window. 
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Figure 7: Region of parameter space in the moMi^ plane where the DAMA signal is consis- 
tent with the constraints from all null experiments. The widest allowed range of dark matter 
masses occurs where 5 = 0. 

The light dark matter window suggested by the DAMA results motivates many new 
directions for model-building and phenomenology with low-mass WIMP candidates. Al- 
though more model-dependent, the implications of light dark matter for indirect detection 
and collider experiments should be explored. It would be interesting to also consider whether 
spin-dependent scattering allows a larger range of dark matter masses to be consistent with 
DAMA and the null experiments. 
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